State reconstruction of a multimode twin beam using photodetection 
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Robust and reliable method for reconstructing quasi-distributions of integrated intensities of twin 
beams generated in spontaneous parametric down- conversion and entangled in photon numbers 
is suggested. It utilizes the first and second photocount moments and minimizes the declination 
from experimental photocount histograms. Qualitatively different forms of quasi-distributions for 
quantum and classical states are suggested in the method. The transition from quantum to classical 
states caused by an increased detection noise is discussed. Momentum criterion for non-classicality 
of twin beams is suggested. 
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I. INTRODUCTION 

Reconstruction of a quantum state from the measured 
data represents one of the most important, and also dif- 
ficult, problems in quantum mechanics. This problem 
has been addressed many times in the past using differ- 
ent approaches. Mathematicians developed the method 
of maximum likelihood (J Q that allows to determine, in 
principle, a quantum state of an arbitrary system that 
suits in the best way to the experimental data accord- 
ing to a given criterion. However, huge amount of ex- 
perimental data needed namely for systems described in 
larger Hilbert spaces (including composite systems being 
in entangled states like states of twin beams) represents a 
serious drawback of this otherwise rigorous method. This 
limits its useful application to only very simple quantum 
systems like those composed of several two-level systems 
or spins. Reconstruction of a quantum state of more com- 
plex systems including those characterized by continuous 
variables inevitably needs different approaches. 

To achieve the preselected level of precision in the 
maximum-likelihood method, the number of needed mea- 
surement repetitions is determined by the smallest values 
contained in a measured distribution. On the other hand, 
usually only large values of the measured distribution 
contribute significantly to physical quantities important 
in the description of the analyzed system. That is why 
a smaller number of measurement repetitions should suf- 
fice. Statistical moments of the distribution represent the 
tool that allows to separate physically important (and 
relevant) information. It is well known that only the 
first and second statistical moments of measured phys- 
ical quantities can be obtained with sufficient precision 
after a reasonable number of measurement repetitions 
On the other hand, the third and higher statistical 
moments can scarcely be exploited in interpreting most 
experiments in quantum physics. 

That is why we have elaborated a method that re- 



lies only on the experimental first and second moments 
of suitable physical quantities. However, these moments 
need not be sufficient for the characterization of a re- 
constructed state, especially in case of entangled quan- 
tum states As the example of twin beams entan- 
gled in photon numbers and analyzed below shows ad- 
ditional parameters are needed for the characterization 
of such states. Values of these parameters can then be 
determined using, e.g., the most common method that 
minimizes square declinations from the experimental his- 
tograms. 



The presented method relies on the measurement of 
photocount statistics. Alternatively, homodyne detec- 
tion based on mixing the analyzed field with a co- 
herent local oscillator may be applied (for review, see 
[5|). However, this method is useful only for fields com- 
posed of one spatiotemporal mode or a small number of 
modes. It allows to reconstruct even the whole phase- 
space quasi-distribution (of field amplitudes) using qual- 
itatively larger amount of experimental data. That is 
why only simpler systems like single-mode fields in Fock 
states with small numbers of photons or squeezed two- 
mode vacuum have been analyzed by homodyne detec- 
tion up to now On the other hand, the method pre- 
sented in this paper allows only for the reconstruction 
of quasi-distributions of integrated intensities, which is 
sufficient for multi-mode fields. The obtained results are 
more reliable as the precision of the proposed method is 
superior above homodyne detection. 



In this paper, the method for reconstructing a state of 
twin beams is presented in Sec. II. Quasi-distributions of 
twin beams are studied in Sec. III. Conclusions are drawn 
in Sec. IV. 
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II. TWIN BEAMS AND THEIR 
RECONSTRUCTION 

The analyzed state of twin beams occurring in the pro- 
cess of spontaneous parametric down-conversion (SPDC) 
is entangled in photon numbers. In an experiment, 
histograms of joint signal-idler (SI) photocount distri- 
butions acquired, e.g., by an iCCD camera are avail- 
able @, H| • Alternatively, time- multiplexed systems with 
avalanche photodiodes semiconductor detector ar- 

rays plj . hybrid photomultipliers p^l [l3j or supercon- 
ducting bolometers pj| can be used to capture photo- 
count histograms. 

The states generated in SPDC are highly nonclassi- 
cal because they are composed of photon pairs. In more 
detail, they are superpositions of states differing in the 
number of photon pairs. Each photon in the signal field is 
accompanied by its twin in the idler field. As the number 
of emitted photon pairs is uncertain, the number of de- 
tected signal photons as well as the number of detected 
idler photons fluctuate. However, because the state is 
entangled in photon numbers, the number of idler pho- 
tojns in an ideal field equals the number of signal photons 
In reality, the analyzed optical field contains both 
photon pairs and single photons originating either from 
pairs in which one twin has not been detected or from 
straylight. Thus the overall SI field is found in a general 
entangled state of a bipartite system composed of two 
infinitely large Hilbert spaces. 

In quantum theory, this state is described by a joint SI 
quasi-distribution P(W S , Wi) of integrated intensities W s 
and Wi (QDII) of the signal and idler fields, respectively 
3, 16]. Alternatively, a normal characteristic function 
Cw defined as 

C w (s s ,Si) = (exp [is s W s + isiWi]}_sr = 

/ dW s dWiP{W s ,Wi)exv[is s W s + i Si Wi] 
Jo Jo 

(1) 

may be used. This is convenient for statistically indepen- 
dent fields as their common normal characteristic func- 
tion Cw factorizes. An appropriate functional form of 
the characteristic function Cw has to be chosen to de- 
scribe properly entanglement in the state. As shown 
below, the use of the first and second moments of in- 
tegrated intensities as parameters is sufficient provided 
that we separate the fully entangled (paired) and noise 
parts of the state. The overall field is then considered 
as composed of three independent components describ- 
ing photon pairs, signal noise photons and idler noise 
photons. Independent components of signal noise pho- 
tons and idler noise photons are assumed to be in the 
usual multi-mode thermal states The corresponding 
normal characteristic functions can be found on the first 
line of Eq. @ below. The form of normal characteristic 
function appropriate to a multi-mode paired field (with 
equally populated pairs of signal and idler modes) orig- 



inating in spontaneous process has been derived in [17j . 
It is written on the second line of Eq. (T2J) below. The nor- 
mal characteristic function Cw of the overall field can be 
expressed as pj], [H[ : 

r i , 1 1 

x • (2) 

(1 - is s B p - isiBp + s s s l B p ) M f ' 

M p gives the number of equally-populated modes with 
the mean number B p of photon pairs per mode. These 
modes of photon pairs form the quantum entangled part 
of the system. On the other hand, there also exist M s 
(Mi) modes with the mean number of B s (Bi) of single 
noise signal (idler) photons. They as noise fields clearly 
belong to the classical part of the system. Whereas five 
independent parameters suffice for the description of a 
classical state with a Gaussian form of QDII, six inde- 
pendent parameters have been introduced for the con- 
sidered entangled state. This increase in the number of 
parameters looks like a consequence of separation of fully 
entangled quantum (paired) and classical parts of the 
beams. In fact, it reflects a completely different struc- 
ture of QDIIs of quantum entangled states compared to 
classical ones with a Gaussian shape (see Sec. III). We 
note that the form of normal characteristic function Cw 
in Eq. ([2]) can be generalized to include also fields emitted 
in stimulated parametric down-conversion [l8j]. 

In the accompanying experiment, photon pairs have 
been generated in non-collinear geometry in type I BBO 
crystal 5 mm long pumped by femtosecond pulses coming 
from the third harmonics of a TkSapphire laser tuned to 
840 nm. Whereas the signal beam has propagated di- 
rectly onto the photocathode of an iCCD camera, the 
idler beam has been reflected on a dielectric mirror (R « 
99 %) first and then impinged on a different area of the 
photocathode. In front of the photocathode, the nearly- 
frequency-degenerate signal and idler photons at 280 nm 
have been filtered (AXfwhm ~ 14 nm). After many rep- 
etitions of the measurement, the camera has provided a 
histogram f(m Sl mi) giving the number of measurements 
containing the specified numbers of detected signal (m s ) 
and idler (mi) photocounts. The number of measure- 
ment repetitions has allowed to reliably determine from 
the measured photocount histogram f(m s ,mi) the first 
and second moments of photocount statistics denoted as 
(m 8 ), (rrn), (m? s ), (mf), and (m 8 rrii). Also the first ((d s ), 
(di)) and second ((d%), (df), (d 8 di)) moments of dark 
count statistics have been experimentally obtained. Sub- 
sequently, the measured moments of integrated intensi- 
ties have been derived according to the relations elimi- 
nating the influence of dark counts: 

(W a ) E = (m a )-(d a ), 
((AW a ) 2 ) E = (m 2 a ) (m a ) 2 (m a ) 

- (d 2 a ) + (d a ) 2 + (d a ), a = s,i, 
(AW a AWi) E = (m a mi) - (m s )(mi) 
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(d s d t ) + (d s }(di); 



(3) 



AW = W-(W) E . 

The signal and idler photons are captured by the pho- 
tocathodc with non-unit detection efficiencies rj s and rji, 
respectively, because of losses in the setup and finite 
quantum detection efficiency of the camera. Once the 
efficiencies r) s and % are known, the moments of inte- 
grated intensities of the fields beyond the crystal can be 
found. They are related to the experimental moments of 
integrated intensities in Eq. ([3]) as follows [3]: 

r)a[(W p ) + (W a }] = (W a ) E , 



rf a [((AW P ) 2 ) + ((AW a ) 2 )] = ((AW.) 



[(W p ) + ((AW P ) 2 )} = 



(AW S AW) E 



VsVi 



(4) 



The last formula in Eq. ((4]) expresses the fact that only 
the field of photon pairs (W p ) creates correlations in the 
signal and idler photon numbers. Integrated intensities 
W s and Wi of the noise signal and idler fields are not 
correlated to the intensity W p of the paired field. Equa- 
tions ((3]) represent five relations for six independent non- 
negative moments (W a ) and ((AW a ) 2 ), a = p,s,i. The 
form of Eq. ((4]) also shows that bounds for possible values 
of these moments exist. 

Possible solutions of Eq. ^ form a parametric sys- 
tem that can be conveniently described by the second 
moment ((AW P ) 2 } of the paired field with the allowed 
values in the range (0,mm[((AW s ) 2 ) /r/ 2 , ((AW,) 2 )/^]). 
According to Eq. Q, the remaining moments are: 



(W P ) = 

(Wa) = 
((AWa) 2 ) = 



(AW s AWi) E 
(W a ) E (AW.AWj) E 



((Awy 2 ), 



Va 

((AWq) 2 )g 

?7a 



+ ((AW^) 2 ), 



VsVi 

-((AW P ) 2 ), a=s,i. (5) 



It can be shown that the solution for the relations (|4]) 
exists only if the following inequality is obeyed: 

(AW s AW,) £ /a - min [((AW S ) 2 ) E , ((AW t ) 2 ) E / a 2 } 



,> 



mm[{W s ) E ,(Wi) E /a] 



(6) 

where a = r)i/ri s denotes the ratio of quantum efficien- 
cies. For typical experimental data, the inequality ([6]) 
puts a strong requirement to the allowed values of quan- 
tum efficiencies rj s and rji. 

Quantum theory of photo-detection Q shows that a 
joint signal-idler photon-number distribution p corre- 
sponding to the normal characteristic function Cw in 
Eq. ([2]) takes the form of a two-fold convolution of three 
Mandel-Rice distributions Q: 



mm n, ,n; 



p(n s ,rii) = ^2 p(n s - n; M S) B s ) 
n=0 

x p(m - n; M h Bi)p(n; M p , B v ), (7) 



where p(n; M, B) = T(n + M)/[n\ T(M)]B n /(l + B) n+M 
and r denotes the gamma-function. Mean photon num- 
bers B a and numbers M a of modes as they were intro- 
duced in Eq. (T2|) are obtained from the moments written 
in Eq. ©: 



B„ 



((AW a ) 
(Wa) 



M (Wg) 2 . 

Ma= ((AWa) 2 ) ' a=P ' S ' Z ' ^ 



A joint signal-idler photocount distribution p c describ- 
ing theoretically the measured photocount histogram / 
is then derived from the photon-number distribution p 
in Eq. ([7]) provided that the detection process is charac- 
terized. An iCCD camera with TV pixels, detection effi- 
ciency r) and dark-count rate D= {d)/N is described by 
the probabilities T(m, n) of having m photocounts out of 
a field with n photons in the form 8] (a = s,i): 



T a (m,n) 



M 



(l-D a ) N »(l-r, a ) n (-iy 



E 



(-1)' 



1 J (1 - Da) 



I Va 

N a 1 - r) a 



(9) 



The formula in Eq. © allows us to express the joint 
signal-idler photocount distribution p c as follows: 

oo 

p c (m s ,m l )= ^ 7 1 s (m s ,n s )T 4 (m i ,n 4 )p(n ;j ,n i ). (10) 

n s ,n,i—0 

The method of least square declinations then mini- 
mizes a function T> defined as: 



V 



\ 



E 



[p c (m s ,mi) - f(m s ,mi)Y 



(11) 



In this minimization, all elements of the experimental 
histogram / are taken into account. 

In the experiment, numbers of detected photocounts 
are monitored in three areas on the photocathode of the 
iCCD camera [§]. The first and second areas are illumi- 
nated by the signal and idler fields, respectively, and typi- 
cally contain several detection events. The last area mon- 
itors noise (straylight and dark counts). The numbers 
of photocounts are collected from typically 10 5 frames; 
each frame arises from detection of the signal and idler 
fields generated from one pump pulse. Subsequently, 
the histogram / giving the number of frames with de- 
fined numbers of signal and idler photocounts is built 
and the first and second moments of photocount num- 
bers are determined. The knowledge of moments of the 
noise monitored in the third strip is used to eliminate the 
effect of this noise to the moments of the measured sig- 
nal and idler photocount numbers. Relations in Eq. ([3]) 
are then applied to derive the moments of experimental 
integrated intensities. In a typical experiment analyzed 
below, the following values of moments have been found: 
(W S ) E - 2.411, (Wi) B = 2.353, ((AW S ) 2 ) E = 0.079, 
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FIG. 1: (a) Function T> giving the declination of theoreti- 
cal photocount distribution and experimental photocount his- 
togram as a function of moment {(AW P ) 2 ). Function T> is 
defined only for the values of ((AW p j 2 ) fulfilling Eq. ®. (b) 
Distribution p sum of the sum n 3 + m of the signal and idler 
photon numbers for {(AW P ) 2 ) = 0.549. 
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FIG. 2: Topo graphs of (a) normally-ordered (s = 1) and 
(b) symmetrically-ordered (s = 0) joint SI quasi-distributions 
P(Ws,Wi) of integrated intensities; s t h = 0.63. Integrated 
intensities W a and Wi are in units of photon numbers. 



((AWi) 2 ) E = 0.095, and (AW s A>V,)b = 0.598. Detec- 
tion efficiencies r/ s = 24.3% and rji = 23.5% have been 
obtained in an independent measurement. 

The dependence of function T> on the 'last free' param- 
eter ((AWp) 2 ) of the investigated state attains a global 
minimum as shown in Fig.JIJa). This minimum is reached 
for ((AWp) 2 ) = 0.549 at the border of the allowed values, 
fn this point, the separation into paired and signal/idler 
single-photon noise fields is such that the numbers of 
modes and their mean photon numbers attain the val- 
ues: M p = 179, B p = 0.055, M s = 8 x 1CT 6 , B s = 320, 
Mi = 8 x 10~ 3 and Bi = 12. This means that the recon- 
structed field contains on average 9.9 photon pairs and 
0.003 (0.1) signal (idler) noise photons. Thus, almost 
99% of the detected photoelectrons have their origin in 
the detection of photon pairs. 



III. QUASI-DISTRIBUTIONS OF INTEGRATED 
INTENSITIES 



According to Eq. (1), QDII P(W s ,Wi) can be deter- 
mined as the inverse Fourier transform of the characteris- 
tic function Cw defined in Eq. @ (s means the ordering 



parameter of field operators [H): 

/>oo />oo 

p 5 (w s ,Wi)= dw' s / dw;p p , 5 (W s - W' s ,Wi-Wl) 

Jo J 

x P a<s (W^Pi, s (W!). (12) 

In Eq. Q12D . QDII P p s of the field of photon pairs takes 
the form |l7|: 



P p ,s(W s ,Wi) 



(WsWi)^'- 1 ^ 2 ( W s + Wi 
■ exp 



ttT(M p )B; 



M„ 



2£L 



x y/—K p , s smc 



Ws-Wi 



S > S t h, 



(w s ws M r- i y 2 

= M f eX P 

r(M p )K p , 5 Dp 
x I Mv -i [2 



Bp. s (W s + W t [ 



DpVWM 



s < s t h! 



(13) 



(1 - s)/2, Dp = ^Bp(B p + 1), K p 



-^p,5 -Bp 

-sB p + (1 - s) 2 /4 and s th = 1 + 2(B p - y/B p (B p + 1)). 
Symbol Im denotes the modified Bessel function and 
sinc(a;) = sin(x)/a;. On the other hand, QDIIs P a s of 
the multi-thermal noise fields occurring in the convolu- 
tion in Eq. (112[) and related to s ordering are derived as 
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1: 

B a s = B a + (1 — s)/2. Provided that the noise multi- 
thermal fields are weak, the QDII P s in Eq. (TT21 inherits 
the features of QDII P p s characterizing the field of pho- 
ton pairs. For the paired field and according to Eq. fp~3|) . 
the QDII Pp^ s maintains its quantum form provided that 
we are close to the normal-ordering of field operators 
(s > s t h). It attains negative values localized in strips 
parallel to the diagonal due to the pairwise character 
of the SI field [see Fig. dja)]. Faithful description of 
these strips with negative values requires the introduc- 
tion of additional independent parameters that complete 
the five parameters of a general Gaussian form. From this 
point of view, the separation of fully entangled quantum 
(paired) part from the rest of the field is extraordinarily 
convenient as it introduces only one additional indepen- 
dent parameter. For lower values of ordering parameter 
s (s < s t h), the superimposed detection noise described 
by the 'ordering of field operators' conceals negative val- 
ues and quantum features (entanglement [l5|) of the SI 
field. This results in a non-negative QDII P ps with a 
smoothed shape [see Fig. Hfb)] that can be successfully 
approximated by a general two-dimensional Gaussian dis- 
tribution function. In this case, the needed five indepen- 
dent parameters can naturally be determined from the 
first and second experimental moments. This is usual for 
any classical field. We note here that a QDII related to 
any value of ordering operator s t h can be derived from 
that related to normal ordering using convolution with 
an appropriate Gaussian function [3j]. However, this pro- 
cedure cannot be inverted for principal reasons Q . That 
is why, we need the QDIIs related to normal ordering for 
full characterization of non-classical states including twin 
beams. 

Considering the general QDII P s in the form of 
Eq. (fT2")l , the threshold value Sth can be obtained from 
the analysis of the first and second moments of integrated 
intensities related to an arbitrary s ordering. The condi- 
tion ([A(W S — Wi)] 2 ) 5th = can be rearranged into the 
formula: 

sth = 1 + 2(0- V7 2 " t) , (15) 

M 8 B S + MjBj + 2MpBp 

M s + Mi + 2M P 
M S B 2 S + MjBf - 2M P B P 
7 ~~ M s + A/, + 2M P 



Inspection of Eq. (15) shows that the threshold value 
s th is below 1 provided that M S B 2 S + MiB 2 < 2M p B p . 
This condition for non-classicality of twin beams can be 
rewritten as 

{(AW S ) 2 ) + ((AWi) 2 ) < 2(W P )._ (16) 
According to Eq. (|16l) . single-photon noise in the signal 
and idler fields has to be sufficiently compensated by the 
paired field to keep non-classicality of twin beams. 

Having a joint SI photon-number distribution p(n s , rij) 
written in Eq. ([7|). the distribution p sum of the sum of 
signal and idler photon numbers with its characteris- 
tic teeth-like structure (only even photon-numbers are 
present) [li| can easily be obtained [see Fig.QJb) above]. 
It provides the simplest experimental evidence of the 
presence of photon pairs in the analyzed field. The pre- 
vailing paired structure of the SI field is also responsi- 
ble for strong sub-shot-noise correlations in the signal 
and idler photon-number difference [{(n s — n,j) 2 ) /((n s ) + 
(m)) = 1 + ((Ws - W t ) 2 )/{{W S ) + (Wi)) = 0.07] [2fJi2. 



IV. CONCLUSIONS 

We have reconstructed a quasi-distribution of inte- 
grated intensities of twin beams using the reliable first 
and second experimental photocount moments and the 
method of least square declinations. Whereas a general 
Gaussian form of the quasi-distribution is suitable for a 
classical field, a more general form is needed for quan- 
tum entangled states for which negative values of the 
quasi-distribution are characteristic. The consideration 
of quantum-classical transition has revealed a momentum 
criterion of non-classicality of twin beams. We consider 
the developed method robust and reliable and as such ap- 
plicable and prospective also in other areas of quantum 
physics. 
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